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RESEARCH HIGHLIGHTS: 

• Equilateral triangular coplanar electrodes are designed and fabricated with a dielectric 
and a protective hydrophobic layer to study the enhanced dynamics and zigzag motion 
of the droplets and compared with results from a square electrode array. 

• The droplet spreading dynamics is found in the underdamped regime. Furthermore, a 
transient force balance approach including the effect of contact line friction is adopted 
to develop expression for the droplet velocity under varying operating conditions and 
successfully compared with the experimental data. 
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The geometry of the coplanar electrodes during electrowetting on dielectric (EWOD) plays a 
crucial role in the dynamics of droplet actuation and motion and the resulting phenomena are 
important in a number of applications. Herein, a strategically placed array of triangular 
electrodes resulted in a zigzag motion of the droplet. Because of the shape and orientation of 
the electrodes, the effective contact line moves in different directions after every switching. 
The displacements, velocities and transfer times are measured to characterize droplet motion. 
A force balance involving the effects of electrocapillary (driving force) and the resistive 
forces, such as the hydrodynamic force and the contact line friction force is developed to 
probe the physics of droplet dynamics on such platforms. The transfer time and the average 
droplet velocities on triangular electrodes are compared with those on conventional square 
electrodes. Further, the spreading dynamics of droplet actuation is probed to elucidate its 
unique characteristics including the enhancement of transport rate for the specific electrode 
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1 . 


INTRODUCTION 


The controlled movement of smáli discrete droplets on a composite hydrophobic- 
dielectric film using the principles of electrowetting on dielectric (EWOD) is being ušed in a 
number of applications leading to the fast-evolving field of Digital microfluidics (DMF). The 
transportation of droplet is governed by an applied electric potential in an array of metal 
electrodes located underneath a dielectric film. There are many applications of DMF e.g., in 
micro-reactors, in biology such as in low volume biological, biochemical and lab-on-chip 
operations like enzyme assay[l], point of care devices[2], DNA based applications[3], 
Chemical reactor bed[4]. Other applications include liquid dispensers[4], electronic hotspot 
cooling[5], electrostatically focused liquid lenses[6], electronic páper display[7], and many 
others[8-10]. 

In generál, the DMF platform is designed in two configurations: closed and open 
platform devices[ 11,12]. In closed DMF devices, droplets are sandwiched between two 
plates. The bottom electrode acts as the activated electrodes and the top one comprises the 
ground electrode. This type of configuration is capable to perform, four different microfluidic 
operations like dispensing, splitting, merging and transporting[10]. On the other hand, single 
plate open platform configuration[11] is suitable for droplet transportation and works better 
in merging & mixing, but lack in performing droplet splitting and dispensing [12]. 

In an open DMF platform, a sessile conductive droplet over the dielectric film covers 
two adjacent coplanar electrodes. Electrostatic actuation is provided by the application of an 
electric potential to one of the electrodes while keeping the other grounded. This leads to a 
reduction in the interfacial energy of the activated electrode, due to the modulation of 
dielectric charges by applied electric field[12]. This reduction in interfacial energy gives rise 
to a change in the contact angle over the electrically activated electrode, while the contact 
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angle remains unchanged over the grounded electrode. This difference between the contact 
angle of the two parts of a droplet may set forth motion of the droplet[13]. From an 
experimental perspective, due to the presence of surface micro defects and surface pinning, it 
requires a minimum actuation voltage (threshold voltage) for droplet actuation[14]. 

Conventionally DMF platforms are designed with square coplanar electrodes for 
droplet actuation[15]. It has been pointed out that DMF performance largely depends on the 
geometry of the coplanar electrodes[16,17]. To achieve smooth transport of droplets with 
enhanced velocity, various research groups háve tried different types of driving electrodes 
like jagged, crescent, twin plate[18,14], etc. Among thcm.intcrdigitatcd jagged coplanar 
electrodes turned out to be the best in different lab-on-chip operations[19]. The jagged 
electrodes perform well for smooth droplet merging and splitting, while crescent electrodes 
report higher droplet velocity, thus increasing the transport frequency for the DMF device. 

In the present study, a new fabrication stratégy is presented by strategically placing an 
array of triangular coplanar electrodes which results in two-directional zigzag motions of the 
droplets with possible benefits in mixing applications. The instantaneous velocities, 
acceleration, and average transfer frequencies over the triangular electrodes are measured. A 
theoretical model based on a force balance approach has been developed taking into account 
the pertinent forces for droplet motion and to analýze the experimental data to elicit important 
information on the role of electrode geometry and effect of contact line dynamics. 

2. EXPERIMENTATION AND METHODS 

The experiments involve actuating the motion of a sessile, conductive liquid droplet 
over the fabricated triangular coplanar electrode array in an open EWOD mode. A single 
array consisting of nine identical electrodes having equilateral triangular geometry with a 
length of each side is equal to 1.5 mm is designed. Each electrode is rotated by an angle of 
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180° with respect to its two neighbors, placed at a uniform gap of 60pm (Fig. 1). The 
triangular electrodes are covered with a thin layer of PDMS as the dielectric and a thin 
coating of Teflon as the protective hydrophobic layer. The details of design and fabrication 
steps are provided in the Supplementary document (Fig. Sl, Fig. S2 of section I and II). 

2.1 EWOD experiments 

The sequential activation of the electrode array for droplet actuation and movement is 
achieved by connecting the fabricated DMF platform to an electronic switching circuit board 
and power supply along with a high-speed camera (Fig. 2). The board consists of Arduino 
Mega with ATMEGA as a microcontroller and is connected to a power supply module 
(Keithley instruments 0-lkV) and the DMF platform through electromagnetic relays for 
sequential activation (Fig. 2(a)). A high-speed camera (IMPEREX: VGA 210 camera) with a 
macro lens (IMPEREX: IPX-VGA210-L) is ušed to capture droplet dynamics at 120 fps (Fig. 
2(b)). 

2.2 The sequence of electrode actuation and droplet movement 

Experiments are perfonned with 4±lpl droplets of lmM potassium chloride solution 
prepared with ultrapure distilled water (Millipore India Pvt. Ltd.). Fig. 3(a) and 3(b) represent 
the schematic top view and photograph of a sequence of droplet motion and the 
corresponding switching logic. The droplet is manually dispensed to ensure that its footprint 
partially overlaps two adjacent triangular electrodes (Fig. 3(a,b)-(i)). A voltage (V) is applied 
to the 2 nd electrode through the activation of relay 1 ensuring positive potential to the 2 nd 
electrode. Once the 2 nd electrode gets activated, the droplet (4pl) starts to move from its 
initial position as shown in Fig. 3(a,b)-(ii) at an actuation voltage of around 175 V. 
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The foot of the droplet after the first actuation rests between the 2 nd and the 3 ld 
electrode. A portion of the droplet may go beyond the 2 nd electrode as shown in Fig. 3(a,b)- 
(iii) as the droplet follows the longest diagonál path of the parallelogram formed by the two 
adjacent triangular electrodes, to reach its finál location (between the 2 nd and 3 rd electrode). 

This is followed by activation of the 3 electrode using the relay 2 such that the droplet starts 
moving towards the 3 electrode. Because of the relative orientation of the two consecutive 
electrodes, the droplet moves in a zigzag path as shown in Fig. 3(a,b) with the arrows 
representing the motion. The video of the droplet motion is provided as supporting 
information and analyzed by an open-source image processing software named ‘Tracker’ at 
four different actuating voltages of 200V, 250V, 300V, and 350V. 

3. RESULTS AND DISCUSSION 

Upon application of an electric field, the droplet starts to spread first, followed by 
movement once a threshold value is reached. The initiation and actuation are addressed 
separately. 

3.1 Droplet spreading dynamics 

The Young-Lippmann equation relates the change in the droplet contact angle upon 

s * s 

application of an electric field as cos 0 V = cos 0 Q H--—-— V 2 where 6 V is the voltage- 

2 -d-r LV 

dependent contact angle, 9 0 is the initial contact angle; s 0 and s r are the relative permittivity 

of free space and that of the dielectric layer having thickness ď, y LV is the liquid-vapor 

surface tension and V the applied voltage. However, the relation provides the equilibrium 
contact angles but does not address the transient States of the droplet during actuation. It is 
mentioned in the literatuře [20] that the transient behavior of a droplet can be categorized into 
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two regimes, námely underdamped and overdamped. In a physical systém, damping is 
produced by processes that dissipate the energy stored in the systém. In the overdamped 
regime, the foot rádius of the droplet monotonically increases and reaches a steady-state, 
while in the underdamped regime the increase in the foot rádius is significantly more and 
decays in a sinusoidal fashion before reaching its equilibrium value[20]. 

The electrocapillary driving force ( F a ) acting at the vicinity of the three-phase contact 
line (TPCL) region of a droplet, F a can be represented by, F = y LV ■ (cos 6 V - cos 9 0 ) = y LV ■ ij 

s * s 

where // is non-dimensional electrowetting number represented as // = —-— —V 2 and y LV is 

2 -d-r LV 

the liquid-vapor surface tension. 

In the overdamped regime, the dominant opposing friction force is the contact line 
friction force[21]. If T od is the characteristic time for spreading of the droplet in the 

r 

overdamped regime, the contact line velocity can be represented as, w cf « ——, where 

T od 

r e = j3-r' is the spreading rádius at equilibrium[20] and r'is the rádius of the equivalent 
sphere having the same volume as that of the droplet. [5 can be derived based on the 
assumption that the droplet retains its spherical cap shape with its volume conserved and 

expressed as/? = — = (1-cos 2 6^) 1/2 •(4-(l-cos6( / )~ 2 (2 + cos6( / )~ 1 ) 1/3 . By balancing the 
r' 

opposing contact line friction force with the electrocapillary driving force, one can write 

C ■ r 

T od =- - , with 2 being the contact line friction coefficient. 

r-v 

In the underdamped regime, the viscous force is negligible as compared to the 
opposing inertial force. By balancing the driving force with the inertial force per unit length, 
it can be shown that the droplet oscillates with a characteristic frequency [22] and can be 
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( V /2 

y ■ v 

represented as co = - 7 . The characteristic time scale to reach maximum 

U-(o 3 J 

deformation can be represented as r ud = n ■ 

The prevalence of overdamped and underdamped regime can be ascertained by 
comparing the two time scales. The ratio of the time scales, denoted by c and defined as 

E, = will indicate the nátuře of the systém during spreading. Any value of f<l will 

^ ud 

represent the underdamped regime and vice versa. The E, values evaluated for the four 
operating voltages (expresed as electrowetting number /;) are presented in table 1. They 
clearly demonstrate that the droplet spreading is well in the underdamped regime with 
dominance of inertial over viscous force for the voltage range ušed herein. 

3.2 Droplet actuation dynamics 

A transient force balance approach involving the actuating force and the resistive forces 
is ušed to evaluate the droplet velocity. The relevant forces are the electrocapillary inspired 
driving force ( F a ) and resistive forces námely, the hydrodynamic shear force ( F/,) and the 
contact line friction force (F c /). The frictional effect due to drag by the filler fluid (air) can be 
neglected. The values of the Capillary number (Ca), Bond number (Bo), and Weber number 
(We) for the droplet of volume 4+1 pl (footprint rádius of the droplet = 0.96mm) at different 
applied voltages are found to be Ca = 1.99x10 5 to 9.50x10 4 ; Bo = 0.257 and We = 5.26x10 5 

o 

to 1.19x10" respectively and the droplet can be considered as part of a sphere. The force 
balance during droplet transport can be written as, 

du „ „ „ 

m- = F a -F h -F d (1) 

dt 
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where, F a = driving electro-capillary force, Z 7 /, = hydrodynamic force, F c i = three-phase 
contact line friction force. 

The two adjacent triangular coplanar electrodes form a parallelogram. Due to the 
change in interfacial tension upon application of a potential difference, the droplet 
experiences a net driving capillary force in a direction perpendicular to the shorter diagonál 
as shown in Fig. 4. The region at the foot of droplet over which this capillary line force acts is 
known as the three-phase contact line (TPCL) region (Fig. 4(a), the are DA of the circle), 
where the droplet, substráte and the surrounding air are in close proximity. 

The capillary line force density acting at the TPCL[15,28]can be represented by, 

fo ~ rLV C OS Oq (1) 

where y LV the surface tension; 9 0 and 9 V are the contact angles in absence and presence of an 

applied potential. Thus, equation (2) can be modified to include the additional effect of the 
electric field as described by the Young’s-Lippmann (Y-L) equation[23] described earlier, 

fa =rLr-(cos^-cos6» 0 ) = ^_Ay 2 (2) 

2 -a 

Since the droplet is moving in a zigzag direction, the TPCL is continually changing. Fig. 4(a- 
c) shows the position of a droplet traveling across two consecutive triangular electrodes 
(respective axes are defined specifically in the schematic for better understanding). 

The x component of the electrocapillary force (f a ) can be represented as 
f ax = f a cos 30 and the y-component as / = f a sin 30°. The net electrocapillary force ( F a ) 
acting along the are DA is therefore, 
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( 4 ) 


K=U(f-) ,+ (fJ dl= f--L n - dl - 

DA 

dl is the unit element of the droplet contour line (are DA of the circle in Fig. 4(a)), n is the 
unit vector normál to the line. Equation (4) can be rewritten as, 


F a = fa ' \ DA n ' dl = YlV ■ ( C0S 6 V “COS 0 Q )'l 


2 -d 


( 5 ) 


Here / is the effective length of the contact line length which is the projection of the 
TPCL over the shorter diagonál of the parallelogram formed by the electrodes. The 
orientation of the triangular electrodes results in changing the direction of the net driving 
force at every switching and results in zigzag motion of the droplet over the electrode array 
and the electro-capillary induced driving force can be written[24,31,32]as, 


F a=l-r L v{ c os6» v -cosé» 0 ) 


l-c„ 




( 6 ) 


Here, 6 V is the instantaneous contact angle and Oo is the contact angle at zero voltage 
(110.5°); with V eff (t) = Vj ] (\-e ,//ír ’)and V eff ° = ^(v 2 -V 2 imn ) ; where V min is the threshold 

voltage for droplet actuation; C eq is the specific capacitance represented by C eq = 

and y LV is the surface tension measured as 70.01 mN.m; 6' 0 is the permittivity of free space, s p 
and s t are the dielectric constants of PDMS and Teflon respectively, d p and d r are the 
thicknesses of the PDMS and theTeflon layer respectively. 

Incorporation of \4//mstě ad of V includes the surface hysteresis effect into the present 
model[27]. The gradual spreading of the droplet increases its rádius of curvature and results 
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in a reduction of the liquid pressure at TPCL, with a concomitant change in the capillary 
force at the TPCL resulting in motion. 

The hydrodynamic friction force can be expressed as, 


K = 


3ju f 7rr(t) 2 u(t) 
4 Ht) 


(7) 


Here, ju f represents the viscosity of the liquid, r is the foot rádius of the droplet, h 

represents the height of the droplet during transportation and u the velocity. The values of r 
and h at different applied voltages are measured and presented in the supporting document 
table Sl. The relation between the contact angle and the footprint rádius of the droplet at 
different applied voltage can be written as[28], 


r(t) = 


3V sin 3 6 V 3 

7r(2-3cos0 v + cos 3 0 V ) 


( 8 ) 


where V is the droplet volume. The height h can be expressed as, 


h(t) = r(t) ■ 


l-cos6^ N 

v sin6 v , 


(9) 


Dynamic wetting takés plače as the droplet starts to move with displacement of the three- 
phase contact line. Such displacement of the contact line leads to a dissipation of energy at 
the molecular level and commonly termed as the contact line friction force[35,36,37], 
denoted here by F c /. The molecular kinetic theory (MKT) expresses the contact line friction 
force as, 


F c , = P£u«) 


( 10 ) 
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where C is the coefficient of contact line friction[26] and P is the perimeter of the droplet. The 


origin of the parameter <T is molecular in nátuře and depends on the nátuře and interactions of 
the molecules at the TPCL region and can be represented as [31], 

r _ k B T 

k°y l 3 (11) 


The evaluation of C is not straightforward and is generally measured in šitu from the 
experimental results. Substituting the expressions of these forces in equation (1) the overall 
goveming equation for droplet motion becomes 


du 

m — = 
dt 


íl-c \ 

eq 1 [V e At)f 


v 


J 


3 J u f 7rr 1 (t)u(t) 
4/7(t) 


PCu(t) 


( 12 ) 


The average velocity ( u avg ) can be represented by numerically solving equation (12) and is 
expressed as. 


u (0 = 

avg x ' 


J u(t)dt 
o_ 


(13) 


Where, t a is the time taken by the droplet to reach steady-state as mentioned in Fig. 5(b) and 
is a function of the applied voltage. Equation (13) can be ušed to evaluate the average droplet 
velocity during electrowetting and compared with the experimentally measured average 
velocity. 

3.3 Evaluation of dynamic characteristics during droplet transportation 

The droplet transportation characteristics are measured at different operating voltages. 
The threshold voltage, at which the droplet start to move is the hysteresis induced minimum 
actuation voltage (V, m „)[30,31]and is found to be 175V for the present experiments. However, 
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smooth translation starts for an applied voltage of around 200V. The visualization of the 
droplet motion clearly indicates that the droplet moves in the direction towards the base of 
the activated electrode, resulting in a zigzag motion (Fig. 3). The advancing edge of the 
actuated droplet is tracked from the extracted sequential images of the recorded video. While 
measuring the displacement of the droplet, the coordinate axis is aligned along the direction 
of motion of the droplet for a single transport thereby taking into account the two- 
dimensional motion over the DMF platform. The changes in the droplet movement 
characteristics for different applied voltage are presented in Fig. 5. Initially, the droplet starts 
to move upon actuation and is marked in Fig. 5(a) as the ‘droplet actuation period’. The 
droplet reaches its peak displacement where no appreciable changes in diameter and velocity 
are observable before the beginning of the next actuation. This is termed as the ‘droplet 
settling period’ in Fig. 5(a). Since the time scales of droplet motion are quite smáli the 
dynamic characteristics are reported in Fig. 5 for durations up to 0.3 s. The charging time of 
the solid-liquid interface (t c ) is found to be very smáli (0.015 s) in comparison to the 
switching time, i.e., t c «t s [32]. Due to the accumulation of a net charge in the liquid-solid 
interface, the resulting electrostatic forces will cause spreading near the contact line, resulting 
in droplet actuation. 

The instantaneous droplet velocities ( u ) are evaluated and presented in Fig. 5(b) at 
different applied voltages with the maximum velocity increasing with applied voltage. During 
droplet transport, the electro-capillary force initially dominates over the frictional forces 
resulting in a peak in droplet velocity. A dynamic equilibrium is eventually reached between 
the electrocapillary force and the resisting frictional forces to reach a region where the 
velocity nearly constant with a very smáli value and is referred to as the ‘droplet settling 
velocity (v,)’. 
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The average transfer frequency (/) is the rate at which a droplet traverses an electrode 
pitch length (e- 1.56mm) during a single transfer[33] (Fig. 5c). Thus,/can be represented 
u 

as / = . The transfer frequency depends on the average velocity of the droplet, which in 

e 

tum depends on the applied voltage. There is a 15% overall improvement in transfer 
frequency over the entire voltage range 200V - 350V. The value of / is indicative of the 
transport performance of the fabricated DMF platform. 

Figuře 6 is a comparison of the experimentally measured values of u avg and that 
predicted fror~ tion (13). The values of //, and<^ are taken to be 8.9xl(T 4 Pa.s and 0.046 

(c) 

Pa.s respectively. The values of č, are obtained by comparing the theoretical (equation 13) 
and the experimental average velocities and are comparable to the reported values in the 
literatuře (0.04-0.08 Pa.s in [31]) highlighting the efficacy of equation (13) for explaining 
droplet transport on a DMF platform. 


3.4 Quantitative estimation of droplet transport on triangular and square electrodes 


The time taken by a droplet to traverse an effective pitch length ( e ) for a specific V is 
termed as the droplet transfer time (t tr ). It is observed that the major difference between the 
conventional square electrode systems and the fabricated triangular electrodes ušed in the 
present study is in the values of t tr for a single transfer. Fig. 7 represents the schematic 
representation of the top view of a sessile droplet over triangular as well as a square electrode 
for a single transfer. The length ( L ) and the inter-electrode gap (p) for both types of 
electrodes are kept the same ( L= 1.5 mm and p =60 pm). The effective length traversed by the 


droplet in a single transfer over triangular activated electrode is defined as, e = 


2 


+ P 


and that for square electrode can be represented as, e =L+ p respectively. Thus the t tr can 
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e 


for triangular and t tr 


e 


for the square electrode. Therefore for 


be expressed as, t u 


identical applied voltage, droplet volume, dielectric materiál, and thickness, t ír should 
approximately scale with e (or e ). As the effective pitch length for triangular electrodes is 
less than the square electrode the estimated increase in a transfer rate over the triangular 
electrode is about 15%. 


The average velocities (u aV g) obtained from the experimental data for triangular 
electrodes are further compared with those from another similar study employing 
conventional square electrodes[34]. It needs to be pointed out that in that study PDMS was 
ušed as the dielectric-hydrophobic layer having a thickness of 15pm and inter-electrode gap 
of 80pm. Thus the comparison will only provide a qualitative indication of the increase in 
transfer efficiency for the present čase. The comparison shows appreciable overall 
improvements in u avg for all the applied voltages (Fig. 8). 

4. CONCLUSION 

The unique design of the fabricated triangular electrodes on the DMF platform resulted 
in a zigzag motion of the droplet due to the continuous changing of the direction of the net 
motive force on the contact line. Significant reductions in transfer time (about 15%) for 
droplet movement along with improvement in average droplet velocities on the fabricated 
DMF platform are observed, quantified and compared to results on conventional square 
electrodes. A model based on a force balance approach between the driving and the resistive 
forces is developed and the resulting equations are ušed to analýze the experimental trends 
and results. The droplet spreading dynamics is thoroughly probed and the values of the ratio 
of the characteristic time scales (4) at the underdamped and overdamped regimes are 
evaluated. The values of č, for the present experimental systém and applied voltage range 
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indicate that the droplet dynamics follows the characteristics of an underdamped regime. The 
values of the contact line friction coefficient are evaluated from the experimental data and are 
found to be consistent with the values reported in the literatuře. These results will be helpful 
in developing micro-actuators with a faster response in lab-on-chip applications. 
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Table(s) 


TABLE FOR MANUSCRIPT 

Augmented Droplet Actuation and Spreading using Triangular Coplanar 

Electrodes 

Table 1: The variation in the values of £, with the electrowetting number, // . 
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FIGURES FOR THE MANUSCRIPT 

Augmented Droplet Actuation and Spreading using Triangular Coplanar 

Electrodes 
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Figuře 1: (a) Schematic view of the designed triangular electrode array with the individual 
contact pads and (b) optical photograph of the fabricated electrode array photomask. 
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Figuře 2: Experimental setup for microdroplet actuation; (a) the Arduino based switching 


logic Circuit, (b) schematic of the setup. 
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Floating electrode 
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I = lst electrode Nc = No connection 

II = 2nd electrode +V = Applied voltage 

III = 3rd electrode Gnd = Ground 


-► Direction of droplet actuation 
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Figuře 3: (a) Schematic and (b) photograph of droplet position at three different switching 
logics during movement; (i) initial droplet position, (ii) position after l st switching, (iii) 
position after 2 nd switching. F a represents the driving force. 
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Top view 



during l st transfer; (a) the initial droplet position (red arrow represents the direction of the 
driving force) (b) position during initiation of l st transfer and (c) position after the l st transfer. 


Positions A, B, and C represent the locations of the advancing region of the droplet over the 
activated electrodes. 
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Figuře 5: (a) the displacement of the droplet at applied voltages of 200V, 250V, 300V and 
350V; (b) corresponding instantaneous velocities after first transport (the lineš are the best- 
fitted curves for the data points); (c) the average transfer frequency at different operating 
voltages (±1-4% error is associated with each experimental value). The line is a guide for the 
readers’ eyes. 


5 






























Figuře 6: Comparison of the predicted and experimental (+1-4% error) velocities at different 
applied voltages. The line is a guide for the readers’ eyes only. 
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(a) (b) 

Figuře 7: Comparison of the distance covered by the droplet over (a) triangular and (b) 
square coplanar electrode array in a single transfer. 
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Figuře 8: Comparison of the average velocities for the fabricated triangular and conventional 
square electrodes in DMF platform. 
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